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f^ , We discuss the impact of the recent KamLAND constraint on the solar anti-neutrino flux on the 

^D I analysis of solar neutrino data in the presence of Majorana neutrino transition magnetic moments 

o ■ and solar magnetic fields. We consider different stationary solar magnetic field models, both 

Oh. regular and random, highlighting the strong enhancement in the anti-neutrino production rates 
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'^ ' that characterize turbulent solar magnetic field models. Moreover, we show that for such magnetic 

L^ ■ fields inside the Sun, one can constrain the intrinsic neutrino magnetic moment down to the level 

j^ . of fi^ < few X 10^^'^fiB irrespective of details of the underlying turbulence model. This limit is 

more stringent than all current experimental sensitivities, and similar to the most stringent bounds 
obtained from stellar cooling. We also comment on the robustness of this limit and show that at 
most it might be weakened by one order of magnitude, under very unlikely circumstances. 
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I. INTRODUCTION 



The KamLAND experiment has recently reached an improved sensitivity on a possible 



electron anti-neutrino component in the solar flux [1|. Their current limit corresponds to 
2.8 X 10^^% of the solar boron Vf^ flux, at the 90% C.L., about 30 times better than the 
previous Super-Kamiokande J2[ and the recent SNO |3[ limits. Solar anti-neutrinos constitute 
a characteristic signature of the spin-flavor precession (SFP) mechanism when non-vanishing 
Maiorana neutrino transition magnetic moments |^] interact with solar magnetic flelds % 

The very flrst evidence of reactor anti-neutrino disappearance published by the Kam- 



LAND collaboration 







has already excluded SFP scenarios as solutions to the solar neu- 
trino problem |8|. However this evidence still leaves considerable room for sub-leading SFP 
effects in solar neutrino physics. It is the latest KamLAND limit on the solar electron anti- 
neutrino flux III], in combination with solar neutrino data, including the recent SNO salt 
phase results p, that ultimately establishes the robustness of the simplest three- neutrino 
oscillation description of the solar neutrino dataJlOJ showing how it is essentially stable 
even in the presence of the SFP mechanism |ll| J48[ . 

Little is known about the detailed structure of solar magnetic flelds and several mod- 
els have been previously used to analyse the SFP conversion. These solar magnetic fleld 
models make different assumptions about the nature (re gula r or random) of solar magnetic 



fields, their magnitude, location and typical scales [12, UJl llJ, ll^ ll^- According to the 



iguiar or ranaom^ 



dynamo mechanism the solar magnetic field is generated close to the bottom of the con- 



vective zone 
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convective zone 



Following this picture, we assume that the field resides within the solar 

3 
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14l |. Moreover, in accordance with the present-day understanding of 
solar magnetic field evolution, the large-scale magnetic field in the solar convective zone is 
followed by a small-scale random component, whose strength is expected to be comparable 
to or even larger than that of the regular one. 

Insofar as the SFP mechanism is concerned, the main difference between random and 
regular magnetic field scenarios is that the former generally give rise to an enhanced rate 
of anti-neutrino production, up to two orders of magnitude when compared to the case of 
regular fields of the same (average) amplitude. This fact has been used in Ref. yj in order 
to obtain more stringent limits on [i^B. Moreover, assuming that random magnetic fields 



are of turbulent origin we were able to extract a limit on the transition magnetic moment 
itself (flu)- We show that, under reasonable assumptions, such bound is comparable to the 
best astrophysical limit that follows from stellar cooling arguments [18[- An alternative 
anaysis of the KamLAND data ^ using the " delta-correlated" model for the solar random 
magnetic field was done in 19|. 

The main purpose of the present paper is to give a more comprehensive description 



of the models and analysis already presented briefly in Ref. [ll|. First to elucidate the 
physics of the enhanced anti-neutrino SFP conversion rates in random field models, compared 
to regular magnetic field scenarios. We show how this follows from the loss of coherence 
of the spin flavour evolution in a fluctuating environment. Second, we show how in a 
solar magnetohydrodynamics (MHD) turbulence model of Kolmogorov type one can use 
the characteristic scaling in this theory in order to obtain a limit on the intrinsic neutrino 
transition magnetic moment fi^. 

The paper is organized as follows. In Sec. lUl we develop a perturbative approach to de- 
scribe neutrino evolution in the presence of convective-zone random magnetic fields, treating 
the magnetic interactions as small correction to the oscillation evolution Hamiltonian. This 
provides a good approximation, fully justified in view of recent KamLAND and solar neu- 
trino data which support the MSW LMA interpretation [10| . We show that neutrinos behave 
as a "Fourier analyzer" reading off only that spectral harmonic of the two-point magnetic 
field correlation function whose space period equals the neutrino oscillation length. Our 
solar magnetic field model is discussed in Sec. IIIH first within the framework of the simplest 
piece-constant correlation cell model with one effective correlation scale Lq, and subsequently 
within a Kolmogorov-type turbulent magnetic field picture. Our discussion illuminates the 
difference between SFP anti-neutrino production rates within random and regular fields, 
as well as the physical meaning of the correlation cell model parameters. The results of 
our neutrino data analysis of the KamLAND u^ limit and future perspectives are given in 
Sec. IIVI We also present in Sec. El a critique of the robustness of this limit and show that at 
most it might be weakened by one order of magnitude, under very unlikely circumstances. 
Finally in Sec. IVII we summarize our results. 



II. NEUTRINO EVOLUTION 



Here we adopt a simplified two-neutrino picture of neutrino evolution, neglecting the 
angle 61^. As we show in the Appendix this is a good approximation, in view of current 
neutrino data. Solar neutrino evolution in the presence of a magnetic field involves then 
only the solai_^niixing angle where 612 = Ogoi = and is described by a four-dimensional 
Hamiltonian 
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where Va = z/^ cos 6*23 — i^,- sin 6^23 (6*23 = Oatm is the atmospheric mixing angle); C2 = cos 20 
and S2 = sin 20; 8 = Am?/4E is assumed to be always positive. Note in this approximation 
the Majorana neutrino transition magnetic moment element fii, = fiea describing transitions 
between neutrino flavour states z/g and Ua coincides with the element fiu characterizing 
transitions between mass eigenstates z/i and z/2; Ve(t) = GFV^{Ng{t) — Nn{t)/2) and K(^) = 
GF\/^{—Nn{t)/2) are the neutrino matter potentials for Uf.^ and UaL in the Sun, given by 
the number densities of the electrons {N(,{t)) and neutrons (A^„(t)). Finally, h± = hx ± ihy 
denote the magnetic field components which are perpendicular to the neutrino trajectory. 

Inside the radiative zone, where the magnetic field is neglected, the evolution of the 
neutrinos reduces to that implied by the LMA MSW oscillation hypothesis. In order to get 
an approximate analytic solution for Eq. (^ in the convective zone it is convenient to work 
in the mass basis. Defining the vectors 



l^L 
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we then express the evolution equation in a block form as 
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(5) 



with \I/ = 1/2 J [Ve{t) + Va(t)]dt. The neutrino eigen-energies Ei^ are defined by 



(6) 



where l^ = (\4 — K)/2 and the minus (plus) sign corresponds to the energy state with i = 1 
(2). The solar neutrino mixing angle in matter is given as 

SS2 



tan 29„ 



(7) 
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Similar expressions for the anti-neutrino eigen-energies and matter mixing angle, Em and 
9m, are easily obtained by changing the sign of the matter potential, V — > —V, in Eqs. (jH)) 
and (0). 

In the region of neutrino oscillation parameters indicated by current neutrino data jlOf , 
^ = <^LMA and 9 = ^lma, the matter effect in the convective zone turns out to be rather small, 
since the ratio of the matter potential V to (5lma is at most around ~ 10"^, at the bottom 
of the convective zone. This allows us to expand the neutrino eigen-energies in powers of 
V/6 



E, 



iL 



t6 



i-y^^ + ^^^ + ^U 



and similarly the neutrino mixing angle in matter may be expressed as 

9^^9 + jS2-^{l + si) tan20 + O ( — 



(8) 
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An analogous estimate gives an upper bound on the mixing angle derivative, \9m\l5 < 
2 X 10"^ and 

" < 5 X 10- r^5:i^V^55^V (10) 



Aiv&± "" " " V ^^u J \ b± J ' 
with similar approximations valid for anti-neutrinos. 

Therefore we can safely neglect all powers of V/6, along with 9m and 9m- In this approx- 
imation the (4 X 4) evolution equation decouples into two (2 x 2) equations 
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describing spin-flavour precession of tlie two mass-eigenstate pairs. In tlie absence of mag- 
netic fields the neutrino eigenstates propagate independently across the convective zone. 
When the magnetic field is switched on one obtains a mixing, characteristic of the spin 

n 

flavour precession mechanism [4|. This turns out to be rather small, a fact that greatly 
simplifies the problem of solving the evolution equation. For this we consider the parameter 
K defined as 

, - ^ - 2 5 X 10- f ^i^V ri!i=^V riiii^l^^V T-^V (13) 

"" ,52 -^■^'^l" VIO-'Vb/ llOOkoJ 1, Am2 ) ^OMeVJ ' '"' 

The smallness of k, allows us to expand the neutrino survival probabilities at the surface of 
the Sun in powers of k, to any given desired accuracy in perturbation theory. The matter 
effect p, |y| is also rather small, and is fully encoded by the phase \Ef, so that when \I' -^ 
one recovers the vacuum result of Ref . |^ . 

From Eqs. (jllj) and ()12|1 we find that to leading order in k, the neutrino mass eigenstate 
probabilities at the surface of the Sun can be written in the form 

WiL\l=Piil-v), W2r\%=PiV, (14) 



W2L\l^=P2il-v), Wm\l^ = P2V, (15) 

(note that unitarity is fulfilled), where the key small parameter t] is given by 

V = ^ [ dh [ dt2Mh)b4h)e-'''^'^-'^^ + c.c.] . (16) 

^ Jo Jo 

Here L is the width of the convective zone and Pi = |z/jL(r = 0.7 Rq)\'^ denote the probabilities 

that solar neutrinos reach the bottom of the convective zone in a given mass state, calculated 

„u.enca.,y . IKe LMA-MSW o^^ p.tu.. ..Ua. lo IKe n,elhod u.d ,n Q. 

The probabilities for neutrinos to be detected in flavour states a = e,fi,e,fi, P^a, are 



then given by 
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Here Pia is the Ui -^ Va conversion probabihty from the surface of the Earth to the detector. 
The phases ^k [k = 1,2) characterize the evolution in vacuum from the Sun to the Earth, 
and are given by 

(21) 



Arr^iD + L^-R^ 



2E 
where D is the Sun-Earth distance, (pk contain the phases due to propagation in the Sun 

up to the bottom of the convective zone and in the Earth. We have checked that (pk can be 

safely neglected for our purposes. 

In summary, neutrino evolution can be understood as follows: 

• neutrinos, generated in the solar core, undergo LMA MSW conversion and enter the 
convective zone as a coherent mixture of UgL and UaL- By numerically solving the 
corresponding (2 x 2) MSW evolution problem we obtain the amplitudes t'eL(O) and 
t'ai(O) at the bottom of the convective zone, which are then used as initial values for 
neutrino propagation across the convective zone. 

• neutrino evolution within the convective zone is considered as approximately vacuum 
oscillations modulated by a small spin- flavour conversion j^ . This is treated in leading 
order in the small expansion parameter k, Eq. p3|). 

• the neutrino survival and conversion probabilities, defined at the solar surface, are 
then evolved to the detector taking into account regeneration in the Earth. 

This way neutrino and anti-neutrino yields are determined and used in the analysis of solar 
neutrino and KamLAND data. 



III. MAGNETIC FIELD MODEL 

The main goal of this section is to calculate the parameter rj in Eq. p6p characterizing 
the solar anti-neutrino production rate to first order in the small expansion parameter k in 



Eq. ()13|) describing neutrino propagation in the framework of different solar magnetic field 
models. 

Different approaches have been used in order to describe neutrino propagation in the 
fluctuating matter density or random magnetic field environments. Since the origin of the 
media fluctuations is unknown, we stick to the simplest "white" noise description, which 
involves just two parameters: the typical correlation scale and its characteristic amplitude. 

The "delta-correlated" model described by 

Vi +r2 



(&±(ri)6x(r2)) = &i 



LoS{ri - ra) 



has been used both for the case of matter density fluctuations 



random magnetic fields 
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21, 


22 



as well as 



2J, |25[ . This choice of correlator has the virtue that it allows 



for analytical solutions to the neutrino evolution equations. It suffers from an important 
limitation, namely that the allowed spatial scales, Lq, of fluctuations should be less than 
typical neutrino oscillation lengths Xosc = ^nE/ Am?. However, it so happens that the 
strongest conversion effect takes place when these two parameters are comparable 26 1. 



Because of this one needs to go at least one step further and consider the so called 
"piece-constant" model. For the case of matter density fluctuations this has been used in 



Refs. 



221 



28l | and for random magnetic fields it was applied in [1^ 
In what follows (Sec. IIII Ajl we adopt this piece- const ant model, whose correlator is de- 
scribed by. 



(&±(ri)&x(r2)) : 



:0 



ri + r2 



if |ri -r2| < Lo 
if |ri -r2| > Lq. 



In other words, the fluctuation correlation function is modeled as a step-function. We refer 
to this as "simplest random magnetic field model" . 

Note that in these models the typical correlation scale is in general a free unknown 
parameter. One may go a step further if one has more information about the nature of 
the assumed fluctuations. For instance, a detailed analysis of how density fluctuations, in 
the form of helioseismic waves, can affect the MSW neutrino oscillations was introduced in 
Ref. [29|. In the discussion we give in Sec. IIII Bl we use the solar MHD turbulence model in 
order to describe the nature of the solar random magnetic fields. These are described by a 



8 



magnetic field correlation tensor 

whose specific form will be given below (Eq. ()27|) ). As a result the separate dependence 
on the correlation scale and the amplitude of fluctuations is replaced by a dependence on a 
specific combination of these parameters, e in Eq. ()35|) . This has the advantage of expressing 
the final neutrino fluxes in terms of a single effective parameter which varies over a relatively 
narrow range. 

A. Simplest random field model 

As already mentioned, current views on solar magnetic field evolution suggest that the 
mean large-scale magnetic field is followed by a comparable small-scale random magnetic 
field component. Such random small-scale magnetic field is not directly traced by sunspots 
or other tracers of solar activity. Dynamically, this field propagates through the convective 
zone and photosphere drastically decreasing in strength. While we lack a direct reliable 
observational estimate of its amplitude, one finds that the ratio of the random to regular 
magnetic field amplitudes may be as large as 50-100, as it is not clear at what stage the 
dynamo mechanism saturates. This issue has been very actively discussed in the literature 
(see, e.g., |l7|,|30| and references therein). 

The simplest random convective-zone solar magnetic field model is obtained by imag- 
ining that the convective zone consists of a set of correlation cells of volume Lq where the 
random magnetic field is assumed uniform, fields in adjacent cells being un correlated. In 
this picture we treat the small-scale random magnetic fields in terms of a single effective 
scale Lq characterizing the size of the correlation cells. We also assume that within each cell 
different magnetic field components transversal to the neutrino trajectory are independent 
random variables with zero mean value (for a more detailed discussion see, for example. 



Ref. [I4). That is, a given realization of the random magnetic field along the neutrino path 
is a stationary random process described by Gaussian statistics. It is also assumed that in 
order to satisfy divergence- less condition, V ■ B = 0, the magnetic field strength changes 
smoothly at the boundaries between adjacent cells. This simplified model seems reasonable 
since we do not expect a strong influence of details of the random magnetic field structure 



near the layers separating adjacent domain cells. 

In order to estimate r] in Eq. (fT^ we divide the range of integration into a set of equal 
intervals of correlation length Lq and average over random magnetic fields in each correlation 
cell, 



N N 



r^kLo pILq 



1 ^' ^' rKM) rlljo 

iv) = 0-^'EE / ^^1 / dt2{b4h)b4t2)e~'''^'^-'^^ + c.c.) , (22) 

^ fc=l l^l J(k-l)Lo J{l-l)Lo 

this results in the final equation 

to>=(i:^^^-0^^^^^. (23) 

where b'^±n is the averaged square of the random magnetic field in the n-th correlation cell. 
We assume that the root mean square (rms) field varies smoothly along the neutrino trajec- 
tory. One therefore clearly sees the cumulative effect characterizing neutrino propagation 
in random magnetic fields, implied by the sum in Eq. (|23|) . For the simple case where all 
rms field amplitudes in different cells are equal to some common magnetic field value the 
above result gets proportional to the number of correlation cells traversed by the neutrino, 
A^ = L/Lq. 

In contrast, for regular magnetic fields the situation is different. The neutrino to anti- 
neutrino conversion probability after traversing the convective zone with a constant regular 
magnetic field of the same amplitude is proportional to 

" - wrm '"^ (^^^^^^) - ^ ^'"'(^ ■ ^) + « ((1^'^) - # • p^' 

that is, what would be expected after passing only one cell. 

Therefore in the random magnetic field case one obtains a sizable enhancement of the 
neutrino conversion probability as compared with the case of a constant magnetic field of 
the same amplitude. This enhancement is explained by the fact that the random nature of 
the magnetic field destroys the coherence in the neutrino evolution. Therefore, instead of 
adding amplitudes one have to add probabilities [31 1. 

Finally, in order to take account of the shape of the rms random field profile we introduce 
a factor 

^—l A-max 

which is S" = 1 for constant rms field and of the order of unity for other sufficiently wide 



spatial profiles, e.g. S ~ 0.579 for "smooth" profile 

10 



121, -S* ~ 0.577 for triangle profile 



Q, 



S ~ 0.782 for Kutvitsky-Solov'ev profile J32|. An alternative estimate of typical shape 
factors comes from the assumption that the magnetic energy density globally follows the 
approximate solar density profile with the scale height H = O.IRq, leading to S* ~ \J H/ L ^ 
1/^3^0.57. 

Using the above definition of the shape factor one can rewrite Eq. ^T^ as 

{n) = ^^^S'^ sin'{6-Lo), (26) 

where ^i^ax i^ ^^^ maximal value of the transverse magnetic field. This shows that different 
random magnetic field models can be effectively characterized by two parameters, b']_^^^S'^ 
and Lq. 

B. Magnetohydrodynamic turbulence models 

A well motivated class of magnetic field models can be considered for which the number 
of parameters can be further reduced by eliminating reference to the effective scale Lq. The 
relevant average ((6+(ti)6_(t2)) of the transverse components of the magnetic fields can 
be characterized by introducing the two-point magnetic field correlation tensor as Mij = 
{b^{rl)b,{r,)). 

For the case of isotropic, homogeneous and non-helical random magnetic fields M^j is 
separated as 

M., = M^ (5., - ^) + M, ^ . (27) 

where r = ri — r2. The longitudinal {Ml) and transverse (M^r) correlation functions depend 

only on the separation distance between the two points, r = |ri — r2|. Given that V ■ B = 0, 

one can write 

MN{r) = ^^{r'ML{r)) . (28) 

2r or 

One can generalize the above definition of the correlation function in Eq. (j?rj) so as to 
cover the case of media which are isotropic and homogeneous only locally. This can be done 
by expressing the correlator in factorized form as J33| 

M^(ri, rs) = F {^-^^^ i^^(|ri - r^l) , (29) 

where the random magnetic field profile factor F depends on the center-of-mass position of 
the two points ri and r^ and the local correlator Ki\i{r) coincides with the one characterizing 
the isotropic and homogeneous case, Mi\f{r), given in Eq. (j2Zj). 

11 



The above form provides a reasonable description for solar magnetic field profiles. Indeed, 
typical solar rms field profiles {F) vary on scales of the order of the density scale height 
H ^ 70000 km. On the other hand K^ involves averaging on much smaller scales 6~^ 
related to the neutrino oscillation length, about hundred kilometers for the LMA-MSW 
case. 

Substituting Eq. ^I7\i into Eq. (jT^ and restricting the coordinates to the neutrino path, 

we get 

iv) = 1^"^ / dzi / dz2MNizi,Z2)cos[2S{zi- Z2)] , (30) 

Jo Jo 

where Zi and Z2 are two points on the neutrino trajectory. 

Substituting Eq. (J29|) into Eq. (J30|l . changing variables and rearranging we get 

iv) = Afi^S^L / d^K^iO cos(25 ■ , (31) 

Jo 

where the integration is extended to infinity because oi 6 ■ H ^ 1 and the shape factor S'^ 

is defined as a continuous analogue of Eq. (J2SI) 

L 

S^ = \ [ i^^^dz . (32) 



^,_1_ jhliz 

^ J b±raax 




In order to estimate Ki\f{$,) we assume that magnetic field evolution in the solar con- 
vective zone is due to the highly developed steady-state MHD turbulence treated within 



the Kolmogorov scaling theory IT], |3^ |35|. In other words, for large magnetic Reynolds 
number, Rm ~ 10^ 17J, the solar MHD turbulence is pumped by the largest eddy motion 
that follows from the interplay between convection and differential rotation of the Sun. The 
size of the largest eddies, Lmax, may be associated with the solar granule size of the or- 
der of 1000 km. Dynamo enhancement subsequently results in the direct cascade of the 
energy of MHD fluctuations to smaller scales J49|. The smallest scale at which turbulent 
motion starts to decay transferring energy into heat is the dissipative scale defined through 

7 r P~'^' ^ '^^ 1 m 

'diss — -'^max-fl'm ^^ ± ill . 

Within the inertial range, /diss < ^ < -^max? the similarity arguments of the Kolmogorov 
theory require the turbulent hydrodynamic (HD) kinetic energy spectrum to scale as E^b ~ 
fc~^/^, where k ~ l/l is the wave number of the eddies of size / [HI, |2J, \^- Similar 
qualitative arguments applied to the MHD case imply that the corresponding Iroshnikov- 
Kraichnan scaling law can be taken as -Emhd ~ k^^^'^ ITJ. However recent theoretical results 

12 



and numerical simulations suggest that the simpler HD Kolmogorov spectrum can actually 
be used even for the MHD case |36l l37| . For our purposes the difference in the power- 
law exponents turns out to be not specially important so that, we take for definiteness the 
Kolmogorov value p = 5/3. Generalization to arbitrary p's is straightforward. See details 
below. 

Let us model the longitudinal correlation function KL{r) as in hydrodynamics 



^^'» - 3 foM li;;:;] ^^'H— '• '''' 



^max 



where T{x) is the gamma-function, K^{x) is the McDonald function of index /i and 6^ is the 
squared rms magnetic field on scale Lmax- The model function correctly reproduces required 
asymptotics of the Kolmogorov theory: 

(i) K,{0) = {bliz)) = 673 , 

(ii) Kl{0) - Kl{1) ~ l^'^ for / < L^ax • 

Using Eq. (|28p and Eq. (j33|) and taking into account that 25 -L^^^ ^ 1 we finally perform 
the integration in Eq. ()3H1 and obtain 50] 

^''^^ 3 r(l/3) (5-L„,,jV3 ^^-^ 52 (5.L^,JV3- 1^-*^ 

which leads, in normalized units, to 

where /xn is the magnetic moment in units of lO^^^yU^, and the ratio e = 
(6/100 kG)/(L^ax/1000 km)i/l 

The ratio e is not known precisely, but one may estimate it assuming equipartition be- 
tween kinetic energy of hydrodynamic fluctuations and the rms magnetic energy at the 
largest (most energetic) scale Lmax [ill 

^— f« — 36) 

Taking f ~ 3 x lO^cm/s and p ~ 1 g/cm J38| we obtain typical amplitude for magnetic 
field fluctuations h ~ 100 kG at the scale L^ax = 1000 km. Assuming that 6 :^ 50 — 100 kG 
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and that the shape factor 5* hes in the range between 0.5 and 1 we estimate that the product 
eS may vary in the interval 0.25 < eS < 1. 

Before concluding this section we note that both types of magnetic field models lead to 
qualitatively similar results for the neutrino conversion parameter {r]) when the magnetic 
field correlation scale is chosen to coincide with neutrino oscillation length Xosc- Indeed, 
in accordance with the Kolmogorov scaling law lP/{6 ■ -Lmax)^^^ — Aosc • &^/-^max = &Aosc ^^ 
the squared rms field at the scale Aosc- This implies that Eq. ()26|) goes into Eq. ()34|) when 
making the replacement (5-L ~ L/Aqsc- This happens because in the context of the turbulent 
magnetic field model neutrinos effectively feel only one scale, namely their oscillation length. 

IV. LIMITS ON NEUTRINO MAGNETIC MOMENTS 

In this section we analyze the limit on electron anti-neutrino flux published by Kam- 
LAND [l|. In Ref. |ll| we showed how this limit makes the determination of neutrino 
oscillation parameters, Am^^; and 9soi, extremely robust against possible existence of spin- 
flavor conversions. This can be used in order to determine the allowed regions of solar 
neutrino oscillation parameters independently of the magnetic field and magnetic moment 
parameters. 

n 

Using the standard x^ procedure (see [8| and references therein) and taking into account 
full set of solar as well as KamLAND reactor neutrino data we have re-determined the 
allowed regions of solar neutrino oscillation parameters, Aml^i and 6soi, within the recent 
version of Standard Solar Model (BP04) 39|. The results are presented as shaded regions 
in Fig. HJ 

Let us first consider the simple random magnetic field model described in Sec. 1111 Al In 
this case neutrino conversion probabilities depend both on the oscillation parameters, Atti^^; 
and 9soh as well as the parameters fJ^lb^^^^S^ and Lq describing the random magnetic field 
model. Using Eqs. (fT^ and ()23|) we have calculated the predicted electron anti-neutrino flux. 
In Fig. ^we show, for a fix value of Lq = 100 km and S"^ = 1, the curves that correspond 
to an electron anti-neutrino yield of 2.8 x 10^^(f)B- It is clear that a better determination 
of the solar mixing angle by future experiments will not substantially improve the limits 
on the parameters /i^&5_ ^"^ ^"^^ -^o which are mainly restricted by the solar anti-neutrino 
flux limit. In contrast note that an improved determination of the solar mass splitting at 
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FIG. 1: Allowed regions for Arn-^^j and Qgoi from solar and KamLAND neutrino experiments 
at 90%, 95%, 99% and 99.73% C.L. The curves are electron anti-neutrino iso-flux contours for 
Lq = 100 km, plotted for different fiiiSb± values indicated in kG, and fixing the anti-neutrino 
production at the current KamLAND limit. For values above 55 kG or so the high LMA region is 
ruled out. 



KamLAND will play an important role in pinning down the magnetic field parameters. 

In order to determine the restrictions on these parameters we have imposed that the 
anti-neutrino yield should not exceed the current experimental bound, 2.8 x IO^'^^b within 
the presently allowed 90% C.L. region of the (Am^^;, 9soi) plane. The results of this analysis 
are shown in Fig |21 The limits on filb^ S'^ versus Lq correspond to different values of 
electron anti-neutrino fluxes. The lower curve represents the current upper bound on the 
product of magnetic moment and magnetic field which follows from the recent KamLAND 
ll. This is compared with the original sensitivity expected by the KamLAND 
401 1 (dashed hne) and with the Super-K jj] bound (dot-dashed line). One can 
see how Fig. |21 quantitatively confirms the expectation that the strongest limit on ^1h\S'^ 
corresponds to the case when the correlation scale Lq is of the same order as the neutrino 
oscillation length, Aosc ~ 100 — 200km. 



bound in Ref. 
collaboration 



As discussed in section IIII B| solar MHD turbulence provides an attractive framework for 
the solar magnetic field model, in which the anti-neutrino production probability depends 
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FIG. 2: Upper bound on \i\\Sh^ for a random field with spatial scale Lq allowed by the 
C. L. region of Am^^^ and Bgoi- The solid line shows the allowed values coming from the current 
KaniLAND limit on 0p^. For comparison we also give the limits that would correspond to ^z,^ = 
0.1% and (pp^ = 0.8 % of (pB- 

only on one extra parameter: /i^e^S*^. As we did above, we first determine the values of 
fi^e'^S'^ that produce an electron anti-neutrino yield of 0.028 % (pB- In Fig. IHlwe have shown 
curves corresponding to different values of finsS. Similarly to the case of the simplest 
random field model one sees that an improved determination of solar neutrino mixing angle 
will not limit fi^^S significantly better than the current constraint. In contrast a better 
determination of the solar mass splitting at KamLAND will be useful. 

Following the same approach as before we have determined the limit on fii,eS taking 
into account the currently 90% CL. allowed region of solar neutrino oscillation parameters. 



As we have already seen in Sec. IIIIB| a reasonable estimate of the allowed range for eS 
is 0.25 < eS" < 1. Therefore, in contrast to the previous models, with random or regular 
fields, we can now, to within a factor of four, extract direct restriction on the intrinsic 
neutrino magnetic transition moment fi,y. This is indicated in Fig. |3J The lowest horizontal 
line represents current limit on the solar electron anti-neutrino flux from the KamLAND [l| 
experiment. On the other hand our bounds on fi,^ are given by the crossings of the lines 
delimiting the dark band with the horizontal line labeled KamLAND. From the Fig. |3]one 
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sees that /ii, < 5 x 10 ^'^jj.b [UI- For comparison the best current laboratory hmit by the 
MUNU experiment {jj^ < 1.0 x IO^^Vb at 90% C.L.) |4l| is also indicated. This should 
be compared with the best astrophysical limit, estimated as fx,^ < 3.0 x lO^^^/ig jl8|. As 
discussed in Sec. |3 under unlikely circumstances this bound might be weakened by at most 
one order of magnitude. This would correspond to the tilted line delimiting the hatched 
band from the right. 
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FIG. 3: Similar as Fig. ^ The curves are iso-flux contours with the anti-neutrino yield fixed at 
the current KamLAND Hmit on (pp^ and corresponding to different values of eSfin. For values 
above 0.15 or so the high LMA region is ruled out. 



An important issue arises here, namely the robustness of the bounds we have obtained 
with respect to different possible choices of the scaling law for the turbulent kinetic spectrum. 
In order to answer this question we have considered the Iroshnikov-Kraichnan model |l7|. 
characterized by the power law p = 3/2 instead oi p = 5/3 that corresponds to the Kol- 
mogorov spectrum. We have found that, in this case, the anti-neutrino yield is higher by 
~ 30%, implying a correspondingly stronger bound on the neutrino magnetic moment. In 
general, values of p lower than 5/3 lead to the same tendency: the smaller p, the stronger 
limit on the neutrino magnetic moment. 

Another question which may be addressed is whether different power-laws could be dis- 
tinguished in future experiments, should solar electron anti-neutrinos ever be measured. 
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FIG. 4: Bounds on fiy for the turbulent magnetic field model described in the text. The horizontal 
lines indicate the bounds on solar electron anti-neutrino fluxes from Super-K and KamLAND. 
The diagonal dark band shows the limit on the intrinsic neutrino magnetic moment, its width 
corresponds to our (Kolmogorov) turbulent magnetic field model uncertainties. The crossing of 
this band with the lowest KamLAND line gives the limit ^^ = 5 x 10~^^ fiB- For comparison the 
vertical line indicates the present MUNU reactor limit. A more conservative and pessimistic limit 
discussed in Sec. EI would give the limit indicated by the hatched band. 

Although the anti-neutrino yield is larger for the Iroshnikov-Kraichnan spectrum than for 
the Kolmogorov one, the expected anti-neutrino spectrum is not significantly different. In 
Fig. the predicted anti-neutrino spectra (normalized to unity) are plotted both for the 
Iroshnikov-Kraichnan and Kolmogorov spectra. For comparison we have also shown the ex- 
pected electron anti-neutrino spectrum for the random magnetic field case with Lq = 100km 
and for the Kutvitsky-Soloviev magnetic field scenario 13|. One can see that there is no 
essential difference between them. Therefore, if a positive anti-neutrino signal is ever de- 
tected, the spectrum shape would not convey definitive information about the turbulent 
energy spectrum. In contrast, different Lq values in random field scenario lead to signifi- 
cantly different energy spectrum predictions. This may help in some cases to distinguish 
these models, but the detailed analysis of this phenomenon is out of the scope of the paper. 
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FIG. 5: Predicted normalized electron anti-neutrino spectrum for the MHD turbulent model with 
Kolmogorov (solid) or Iroshnikov-Kraichnan (dashed) scaling exponents and the simplest random 
magnetic field model. For comparison we also indicate the expectation for the a regular magnetic 
field model. 



Note that in the framework of our turbulent picture there is no time dependence in 
the magnetic field correlators, so that the resulting solar neutrino fluxes will have no time 
variation. This is a reasonable aproximation since the neutrino data we use are taken 
over periods larger than characteristic turbulent time fluctuations of solar magnetic fields. 
Possible neutrino flux variations in the Super-Kamiokande data have been considered 42, 
43|. We note however that in the framework of the spin flavour precession of Majorana 
neutrinos, irrespective of the magnetic field model adopted, the expected variations in the 
electron and muon neutrino fluxes should not exceed the current electron anti-neutrino flux 
bound 2.8- 1O~^0B- Therefore time variations are expected to be too small, given the current 
statistical errors in the solar neutrino data. 



V. A CRITIQUE OF THE ROBUSTNESS OF THE BOUND OBTAINED 

We now discuss in more detail the possible weaknesses of the method we have employed 
to constrain the neutrino magnetic moment in the case of turbulent magnetic fields. The 
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bound is given from Eq. (jM|) and depends on Lmax- We have taken Lmax ~ 1,000 km, 
arguing that this is the size of the granules on the surface of the Sun. One might argue, 
however, that the scale L^ax niay be associated with the size of the convection cells which 
is about a few xlO^ km. In such a worse-case scenario the neutrino conversion would be 
smaller, thus weakening the limit. However at most this would weaken our bound by a factor 
(few X 10)^/'^ ~ 3 due to the mild power dependence 1/3. Moreover, the size of the largest 
turbulent eddies, Lmax, is still an open question, see for example Ref. |l7|. 

The bound on the neutrino magnetic moment also depends on the shape factor 0.5 < 
S < 1. Such shape factor S accounts for the fact that the random magnetic field may not 
be located in the whole convective zone (see discussion after Eq. (j25p ). The bound on the 
neutrino magnetic moment assuming that the shape factor is S=0.5 is a conservative bound 
and corresponds to the assumption that only 1/4 part of the convective zone is effectively 
filled by the random magnetic field {S"^ = 0.25). 

The strength of the random (as well as the regular) magnetic field in the solar convective 
zone is unknown. There are different theoretical estimates of its value based on available 
experimental data. In order to derive our conservative limit on the neutrino magnetic 
moment we have used the value 50kG. A smaller value, say 20kG, for the strength of the 
large-scale regular magnetic field in the convective zone would result in a limit only 2.5 
times weaker than the value we obtained. However we stress that the small-scale random 
magnetic field is expected to be comparable or even larger than the regular one. The value of 
50 kG for the random magnetic field does not contradict the helioseismic data, in particular 
the analyses of SOHO MDI or GONG data. For example, the observed asphericity of the 
inferred sound speed may be exp lained by a magnetic field at the bottom of the convective 



zone of strength 200-300kG 



45[. We mention also that recent results on the modelling of 



the sunspot generation suggest that the strength of the toroidal field at the base of the solar 
convective zone must be on the order of 100 kG or so 4^ in which case the bound would 
get stronger. 

All in all, taking all of the uncertainties mentioned above into acount we find that, in the 
worst case, they will be able to weaken our bound at most by one order of magnitude (see 
Fig. 13)). Still the resulting bound is better than the current direct experimental bound. This 
situation would, in our opinion, be extremely unlikely as it would require a conspiracy. 
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VI. CONCLUSIONS 

We have considered neutrino spin-flavor precession in stationary solar magnetic fields, 
treated perturbatively with respect to LMA-MSW neutrino oscillations. We have discussed 
the impact of the recent KamLAND constraint on the solar anti-neutrino flux on the solutions 
of the solar neutrino problem in the presence of Majorana neutrino transition magnetic 
moments. This leads to strong limits on neutrino spin-flavor precession, involving fiyB. We 
analysed these constraints for a number of models of solar magnetic flelds, both regular as 
well as random. We found that for a turbulent solar random magnetic fleld model, one 
can flnd a rather stringent constraint on the intrinsic neutrino magnetic moment down 
to the level of fii, < few x lO'^'^fi^, similar to bounds obtained from star cooling. Such 
magnetic moments would have no effect in the detection process, given current experimental 
sensitivities. We also discussed how, in the worst possible case, these limits might deteriorate 
at most by one order of magnitude. We note also that for the complementary case where 
there is no spin flavor precession in the Sun and the only effect of the neutrino magnetic 
moment happens in the detection process |44| the current sensitivity is weaker than found 
here J47[. Therefore we conclude that turbulent solar magnetic flelds provide an enhanced 
sensitivity to very small neutrino transition magnetic moments. We have also shown how our 
result is rather insensitive to the details of the assumed model of turbulent random magnetic 
fleld. We verifled this explicitly for Kolmogorov's and Iroshnikov-Kraichnan spectra. Should 
solar anti-neutrinos ever be detected, it is unlikely that this would be of help in testing the 
intrinsic scaling law B(L) oc L^^^ characterizing turbulence. 
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APPENDIX A: 3-i/ DESCRIPTION OF D^ PRODUCTION IN SOLAR RANDOM 
MAGNETIC FIELDS 

To justify the validity of the 2-z/ scenario one can generahze it to the 3-z/ case. The prob- 
abihty of the solar electron anti-neutrino appearence at the surface of the Earth (averaged 
over the Sun-Earth distance) takes the form 

Pee ~ |/X12|'/i2{c?3C?2P2l(0.7)+C?3S?2Pil(0.7)} 

+ |/i23r/23{c?34^3L(0.7) + 4P2L(0.7)} (Al) 

+ |/ii3|2/i3{44P3i(0.7) + 4Piz.(0.7)} . 

Here iXjk are the transition magnetic moments; Sjk and Cjk are corresponding sin and cos of 
the 3-z/ mixing angles (^12 - solar, 6'23 - atmospheric, ^13 - reactor mixing angles); Pjiif = 
0.7Rq) are the solar neutrino probabilities at the bottom of the convective zone. The 
integrals Ijk are the generalizations of Eq. (|TH|l 

Ijk = l [ dh [ dt2(6+(ti)6-(t2)e-'''^^(*^-*^)+c.c.), (A2) 

^ Jo Jo 

where Sjk = Am|^/4£', Am'jf^ = ml — rrij. To derive the above equation we have used the 

perturbative approach, as in SecjUl leaving only terms quadratic in magnetic moments in 

the final probabilities. 

We can notice fest that, at a ve.y good app.oxnnatio,, one has ^,(0.7) = .?3 « 1 
(see |4^ and references therein). The u^^ yield in the solar core propagates as in vacuum 
because matter effects are strongly suppressed, as ^23 ~ ^13 ^ Vmsw- 

Therefore the second and third terms in Eq. ()A1|) are proportional to sin^6'i3. In order 

to estimate Ijk we consider the simple random magnetic field model discussed in Sec illl Al 

^2^2 ^ 
^jk = -?2— l-sin^('^ifc-^o) (A3) 

and the turbulent magnetic field model discussed in Sec. IIIIBI 

h2q2r r 

yjjk^vaaxj 

Taking into account that solar and atmospheric mass splittings have significantly different 
scales, Am^2 ^ ^"^23 ~ ^^^13) '^6 obtain in both cases (random and turbulent) that 
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-^23 ~ -^13 ^ lu , which means that neutrino spin-flavor conversion in the channels 1 — 3 and 
2 — 3 is strongly suppressed with respect to the 1 — 2 channel. 

We may conclude that possible constraints on |/i23p and |/Ui3p are very poor. They are 
strongly suppressed by two facts: 

• 5^3 < 0.054 (3(T C.L.) [lo|, 

• 0.018 < Amj^/Amj^ ^ AmfJAml^ < 0.053 {3a C.L.) [lo|. 

Therefore the contribution of other channels involving fi23 and /ii3 to electron anti- 
neutrino production is strongly suppressed both directly by the small value of the angle 
6*13 and by the small ratio of solar to atmospheric squared mass differences A^^^/A^j^. As 
a result we adopt the 2-z/ picture, characterized by a single component of the transition 
magnetic moment matrix (/ii2) as a very good approximate description of anti-neutrino 
production. 
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There is a loophole to this argument, namely the case of Dirac neutrinos. This is certainly 
a particular case of the SFP Hamiltonian 141] for which the anti-neutrino argument does not 
apply. However for this case the gauge theoretic expectations for the attainable magnitudes 
of the transition moment are significantly lower than those for Majorana neutrinos. 
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